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We consider a spin- 1/2 ladder with a lerromagnetic rung coupling J± and inequivalent chains. 
This model is obtained by a twist (9) deformation of the ladder and interpolates between the 
isotropic ladder (9 = 0) and the SU(2) ferromagnetic Kondo necklace model (9 = tt). We show 
that the ground state in the (9, J±) plane has a finite string order parameter characterising the 
Haldane phase. Twisting the chain introduces a new energy scale, which we interpret in terms of a 
Suhl-Nakamura interaction. As a consequence we observe a crossover in the scaling of the spin gap 
at weak coupling from A/Jn oc J±/Jn for 9 < 9 C ~ 87r/9 to A/Jn oc (Jx/Jy) 2 for 9 > 9 C . Those 
results are obtained on the basis of large scale Quantum Monte Carlo calculations. 

PACS numbers: 75.10.Pq, 71.10.Fd, 73.22.Gk 



Low-dimensional quantum magnets are fascinating ob- 
jects from both experimental and theoretical points of 
view. Spin- 1/2 ladders have been widely studied and 
interpolate between the physics of one-dimensional an- 
tiferromagnetic (AF) spin chains and two-dimensional 
systems In the one-dimensional (ID) case, there 

is an important mapping between spin-1/2 Heisenberg 
AF chains and Luttingcr liquids [2| which allows to 
treat such chains by means of exact fcrmionization and 
bosonization methods, resulting in a well-understood 
gapless phase Jj|. Coupling identical chains to form a 
spin ladder is however not a trivial task from a theoreti- 
cal point of view [1,01 . Indeed, the coupling is a relevant 
perturbation and, up to logarithmic corrections, opens a 
gap proportional to the interchain coupling J± @, Q- 
In this paper, we will focus on the opening of the spin 
gap for the case of two inequivalent chains coupled with 
a ferromagnetic rung coupling J±_ < 0. This model is 
dubbed the Spiral Staircase Heisenberg Ladder: 


























\ 



a) 



b) 



c) 



H = J,, J2 (Si.< • Si )i+ i + cos 2 (9/2) S 2 ,, ■ S 2 , m 
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Here S a ^ is a spin-1/2 operator on leg a and lattice site i. 
J\\ > sets the energy scale and the interchain coupling 
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FIG. 1: (Color online) (a) Sketch of Spiral Staircase Heisen- 
berg Ladder, (b) View of the model from the top. (c) For 
9 — 7T, the model maps to the ID SU(2) ferromagnetic Kondo 
necklace model M. 



is taken to be ferromagnetic J± < 0. Geometrically, this 
model may be interpreted as a result of twist deforma- 
tion of a 2-leg ladder (Fig.[T^,) with twist performed along 
one of the legs. Such a spiral structure is characterized 
by the angle 9 (see Fig. [Ud) and interpolates between 
the isotropic ladder (6 = 0) and a ferromagnetic SU(2) 
Kondo Necklace 0] model (9 = tt) 0, M, M, El ■ A moti- 
vation to study this specific geometry comes from the fact 
that a realization of the model schematically presented 
in Fig. [TJ: was synthesized as a stable organic biradical 
crystal PNNNO [Hj] . Possible candidates for realizations 
with twist angle < 9 < tt might be found in the families 
of molecular chains decorated by magnetic radicals. 

In the strong coupling limit, |Jj_/J||| 3> 1, the model 
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maps onto the spin-1 Heisenberg chain with effective ex- 
change interaction J e g = Q (l + cos 2 (0/2)). This phase 
has a spin gap given by A H /J e & = 0.41048(6) [HI 
and is characterized by a non-local string order parame- 
ter [l6| (see a recent discussion in fj~7j): 



(a) 



(O s (n)) = (^ n exp 



n +n 



(2) 



with Sj 



Si j -\- 5*2 j • The expectation value picks up the 
hidden antifcrromagnetic ordering. At weak couplings, 
the analysis depends on the twist angle 9. For small 
twist angles (i.e. close to the isotropic case), one can rely 
on the bosonization and numerical results of Refs. 0, @ 
which yield a spin gap proportional to \J±\ up to loga- 
rithmic corrections. On the other hand, at 9 = tt the 
spin velocity on the second leg vanishes thus inhibiting 
the very starting point of Ref. [y. Alternative approaches 
such as a mean-field theory based on a Jordan Wigncr 
transformation, which yields the correct result for the 
isotropic ladder, predicts a spin gap A oc J\j Jy at 9 = tt 
[l8| . A flow equation calculation has recently been car- 
ried out for the SU(2) Kondo necklace model [lj|, (i.e. 
6 = tt in Eq. (TT])) and is interpreted in terms of the 
onset of a spin gap irrespective of the value of Jj_/J\\- 
To disentangle this situation, we have performed large 
scale quantum Monte Carlo (QMC) simulations of the 
ferromagnetic spiral staircase model. Two variants of 
the loop algorithm [2(| were applied. For the string or- 
der parameter and the spin-spin correlation functions, we 
used a discrete time algorithm and extract the spectral 
functions via stochastic analytical continuation schemes 
[2ll . [22l | . For the spin gap calculation, a continuous time 
loop algorithm was used, where the gap is calculate d by 
a second moment estimator of the correlation length [l5| • 

Our results for the spin gap in units of J c ff in the (6, Jj_) 
plane are plotted in Fig. [2] Enhancing the twist angle 
from 9 = to 9 = tt/2 leaves the spin gap, measured 
in units of J e ff, next to invariant thereby showing that a 
small twist is an irrelevant perturbation [23j . For larger 
values of 9, A is suppressed, and in the limit 9 = tt the 
approach to the Haldane value in the limit J± — > — 00 is 
surprisingly slow. At small values of |J_l/Jii|, and 6 = 
we reproduce the results of Ref. namely A oc J± (see 
Fig. [2}}). Here and in what follows, we neglect logarith- 
mic corrections in our discussion. Fig. [Jb shows that this 
weak coupling behavior of the spin gap is sustained up to 
9 < 9 C ~ 87r/9. Beyond this critical angle (24|, the data 
allows for different interpretations. Let us concentrate 
on the twist angles 9 = 8tt/9 and 9 = tt. A linear ex- 
trapolation of the data would lead to the vanishing of the 
spin gap at a finite critical value of J±. However, in this 
parameter range, we find a finite string order parameter 
(see below), incompatible with a gapless phase. As sug- 
gested by a Jordan- Wigner mean-field analysis [IH, we 
instead assume the existence of an inflection point and 
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FIG. 2: (Color online) (a) Spin gap A(Jj_) as a function of 
|J±/J||| for different twist angles 9. The gap is rescaled by 

Jcft = (l + cos 2 (6'/2)) such that in the large- 1 J± |- limit, it 
converges asymptotically toward the Haldane gap of a spin-1 
chain. At weak couplings, we have carried out QMC simula- 
tions up to /3J|| = 2500 and 2 x 512 spins to ensure size and 
temperature convergence. Inset: zoom on the weak coupling 
region, (b) Results for spin gap on a semi-logarithmic scale. 



fit the data to a quadratic form in the limit Jj_ — > (see 
inset of Fig. [2k). Let us note, however, that we cannot 
exclude the possibility of an exponential scaling. 

The scaling of the spin gap at 9 > 9 C implies a rapid in- 
crease of the spin correlation length £ cx </y/A. For 6 = it 
and J±/J\\ = —0.5, spin correlations decay exponentially 
with characteristic length scale £ ~ 115 (see Fig. [3]). At 
J± j J|| = —0.3 no sign of exponential decrease is apparent 
on the considered 2 x 800 lattice. This is consistent with 
a spin gap decreasing as J±/J\\ (or quicker). Indeed, such 
as scaling leads to £ > 300 which is comparable to the 
largest distance L/2 = 400 accessible in our simulation 
of a 2 x 800 lattice. 

On length scales \i — j\ < £ the spin-spin correlation 
functions follow a slow power law. In particular the data 
of Fig.[3]at J±/J\\ = —0.3 are consistent with S(\i — j\) oc 
(— — jl -1 / 3 . At 9 = 7r, the effective interaction 
on the second leg is set by the Suhl-Nakamura (SN) [25[ 
interaction (26[. In second order perturbation theory, 
without attempting any self-consistent calculation, this 
interaction takes the form in Jsn(q) oc Jj_Xs(q, ^> = 0) 
in Fourier space. Here, Xsili w = 0) is the spin sus- 
ceptibility of the spin 1/2-chain. A first step towards 
a self-consistent treatment is to allow for a gap, A, in 



3 



■JJJ, = -0.3 
J±/J« = -0.4 

j ± /j, = -o.5 m 

J±/J t = -0.6 




(b) 





i 


1 ; 








rJ ± /J, 


= -o.3 m 






= -0.4 u X J 






= -0.5 mi 




■ j±/A 


= -0.6 IB-: 
I 


1 . 



FIG. 3: (Color online) Spin-spin correlation function on both 
legs for the Kondo necklace model (6 = 7r) at different cou- 
plings J±/J\\ on a 2 x 800 lattice. Simulations are carried 
out at /3J|, = 7000 (J±/J\\ = -0.3,-0.4), /3J {] = 5000 
{J±/J\\ = -0.5) and /3J|| = 2000 (J±/J|| = -0.6). 
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FIG. 4: (Color online) Dynamical spin-spin correlations at 
J±/J\\ = 1 for the ladder system (9 — 0) and the Kondo 
necklace model (9 = n). Here we consider a bonding combi- 
nation of the spins across the rungs. (J3J\\ = 200, L = 100) 



Xs{q,^ = 0). Thereby and in real space we expect SN 
interaction to have a range set by £. We interpret the 
above mentioned very slow decay of the spin-spin corre- 
lations on both legs and on a length scale set by £ as a 
consequence of the SN interaction. The SN interaction 
at 9 = 7T sets a new low-energy scale in the problem, cor- 
responding to the slow dynamics of the spins degrees of 
freedom on the second leg. Due to the ferromagnetic cou- 
pling between the chains, this slow dynamics will equally 
dominate the low energy physics of the spins on the first 
chain. This new energy scale is also apparent in the dy- 
namical spin structure factor S(q, oS) plotted in Fig. |4] 
As apparent, a narrow magnon band emerges as the an- 
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FIG. 5: (Color online) (a) String order parameter O s and Oh 
as a function coupling J±/J\\ and several twist angles. For 9 = 
8tt/9,tv finite size effects are still present for the considered 
L — 800 lattice in the parameter range |J_l/Jii| < 0.5. For 
> 1-0 the system size L — 400 is sufficiently large 
enough to guarantee convergence. Simulations are carried 
out up to f3J\\ = 7000. (b) Finite size scaling of the order 
parameters for the parameter sets J±/J\\ = —0.2, 9 = 87r/9 
(blue) and J±/J\\ = —0.3, 9 = n (red). The data for Oh are 
fitted to the form: Oh a L _t * exp( — L/£). 



gle 9 grows from to it. To lend support to the inter- 
pretation in terms of the SN interaction, we have checked 
with exact diagonalization methods that the width of the 
magnon band at 9 — it indeed scales as J±/J\\ in the 
weak inter leg coupling limit (data not shown). In the 
vicinity of 9 = ix, we hence expect that the low energy 
effective model is given by a spin-1 Hcisenberg chain with 
exchange coupling set by the SN interaction. Assuming 
the validity of this low energy model, we predict a spin 
gap which scales as Js n oc J\ / Ju . 

The above arguments and data suggest that irrespec- 
tive of the twist angle and coupling Jj_, the ground state 
of the model corresponds to a Haldanc phase. 

We confirm this point of view by computing the string 
order parameter O s = {O s (n))\ n=L / 2 on a 2 x 800 lattice 
(see Fig. [5K), which is finite in the Haldane phase 
Strictly speaking, this is not a sufficient condition to as- 
certain the Haldane physics since we also need to show 
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thermodynamic limit (when both O s > and Oh > 0, 
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an Ising order is present [l6[). 

In the region where the correlation length £ exceeds the 
lattice length, finite-size effects are present (see caption 
of Fig. [5]). In particular when the lattice size is smaller 
than the correlation length, both Oh and O s take non- 
zero values, since the very slow decay of the spin cor- 
relations mimics Ising type order. As the system size 
grows beyond the correlation length, Oh decreases expo- 
nentially whereas O s in enhanced. Those size effects are 
explicitly shown in Fig. [5b at Jj_/J\\ = —0.2, 9 = 8it/9 
where L>^ and J±/Ju = —0.3, 9 = tt where our maxi- 
mal system size barely exceeds the estimated correlation 
length. Taking those size effects into account, we con- 
clude that in the thermodynamic limit, only the string 
order parameter O s is finite in the whole (9, Jj_) plane. 

In conclusion we have established that the ferromag- 
netic spiral staircase is a Haldanc system, irrespective on 
the twist 9 and coupling constant Jj_. In the weak cou- 
pling region, twisting the ladder introduces a new low 
energy scale which we interpret in terms of a SN interac- 
tion. As a consequence and for 9 > 9 C ~ 87r/9, we have 
provided numerical data showing that at weak coupling, 
the spin gap decreases quicker than the linear Jj_ behav- 



ior of the 2-leg ladder (9 — 0). Analysis of the data is 
consistent with the picture that, for 9 > 9 C , the spin gap 
tracks the SN scale and is hence proportional to Jj_/J\\- 
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